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1 Introduction 



The properties of positronium provide some of the most precise tests of QED. Both the 
experimental and theoretical considerations have reached a very high level of precision, 
requiring for the latter the computations of many-loop diagrams 1.. In the present talk, 
we will address one particular aspect of the current QED bound state models, namely 
the factorization between the bound state dynamics and the annihilation process. The 
present study is motivated by the recurrent contradiction between factorized models and 
Low's theorem^. Our central result is a simple alternative method that allows an exact 
non-perturbative treatment of binding energy (BE) effects, and which produces bound 
state decay amplitudes with correct analytical behaviors. 

Having in hand a formalism in which the BE is treated non-perturbatively is especially 
interesting in the context of quarkonium physics. Compared to positronium, quarkonium 
binding energies are much larger and in addition, they cannot be related to the coupling 
constant. Implications for a number of quarkonium puzzles are presented. 

The present talk is based on a series of papers to which we refer for details and 
references. 

1.1 Basic formalism: factorization 



Positronium 




Figure 1 : Basic factorization of bound state decay amplitudes. 

Historically, the first theoretical model designed in the forties to compute decay rates of 
para and orthopositronium ( J = and 1) states was (see Fig.l) 

r {Ps (J) ^ n-f) = IM"^ {^Vrel(Tscatt. (e+e" ^ ^7))^^^^^^ (1) 

which amount to replace the initial flux factor of the scattering cross section by the prob- 
ability of e^e~ contact inside the bound state, i.e. the Schrodinger wavefunction at zero 
separation \4>o\^ = a^m-^/Syr with m the electron mass. Using this naive factorization 
formula, it was computed 

OcTTl __ 

r [p-Ps — > 77) = — ^ — Pirenne and Wheeler 0] (2) 
T {o-Ps 777) = ma^— Ore-PowellH] 

Currently, a more developed type of factorization is used: the bound state decay 
amplitude is constructed as a three-dimensional convolution integral of the momentum- 
space wavefunction with the amplitude for the free constituents to scatter into the final 
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state (M is the positronium mass, 7 the BE): 

M (Ps (J) ^ n7) = V2M I ^^^scatt. {e^,^e\^' - -7)^^^,^^ (3) 
with (j) (k) = (pOTT^^^^^, 72 = - ^ m'^a^/A 

Using this basis, theoretical computations have reached a great level of precision: 

Tp.ps = ^{l + 5Tp.p,), r^p, =ma^^i^^^— ^(l + 5ro_p,) (4) 
with the perturbative series (see and references in [3]) 



bVp.ps = -Ap- + 2a^ In - + Bp^ - — In^ - + Cp— In - + 5a^^ 5) 

, a 1 3a'^ , n 1 ^ , 1 . 

5ro_p, = -^o 5- In - + - — In^ - + a— In - + ^5^^ (6) 

vr J a vr^ 2vr a vr a vr^ 

with 



Ap = b- vrV4 2.5326 Aa = 10.286606 (10) 



547 = 0.274(1) 



Sp = 5.14(30) i?o = 44.52 (26) /"_niqn\ 

Cp = -7.919(1) Co = 5.517(1) 057-u.iyu; 

1.2 Analyticity? Why and how 

From a quantum field theory perspective, the positronium has the quantum numbers of 
a neutral boson. As such, its radiation emissions are constrained to obey very general 
requirements. From QED gauge invariance and quantum field theory analyticity, the low- 
energy end of the photon energy spectrum is unambiguously predicted. This is the content 
of Low's theorem Note that the analyticity requirement invoked here has nothing to do 
with non-analytical terms (In a) in the series expansions © and © , but instead refers to 
a general property of decay amplitudes as functions of external photon energies, following 
from microcausalityinj. 

To state Low's theorem precisely, let us take the amplitude for orthopositronium to 
three photons 

A^(o-Ps^777) =/(^7-") (8) 
Gauge invariance and analyticity imply that the Laurent expansion of the amplitude for 
small photon energy has no singular and no constant term, and thus vanish linearly: 

M {o-Ps 777) = O (w^) near = (9) 

As we now show, factorization-type approaches are in contradiction with this require- 
ment. Using the naive factorization model Ore and Powell computed the energy 
spectrum [S] 

(iP , „ , 2ma^ ^ , s , 

— (o-Ps^777) = -^Si(xi) (10) 



X {2-xf 



(1 — x) (1 — X 



,2 



(2 - xf 



In(l-x) (11) 
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with xi = 2uj^/M the reduced photon energy. For very soft photon, this spectrum behaves 
as 

5 

n[x) = -X + (x^) near x = (12) 
3 

while an O (x^) behavior is required from Q. The problem originates in bremss-trahlung 
radiations contained in agcatt. of since its corresponding amplitude behaves as 

Mscatt. (e+e- -^nj)=0 {uj-') + O (c^O) + O K) + ... (13) 

While the O (w^^) term is cancelled in the limit Vrei — *■ (selection rules), the constant 
one is not disposed of. Since it is a bremsstrahlung radiation, it is as unphysical as an IR 
divergence for positronium. 

The convolution- type model © does not solve the problem, as can be shown using the 
language of dispersion relations. Starting from a four-dimensional loop model, and taking 
into account only the vertical cut (see Fig. 2) in the imaginary part, the dispersion integral 
is precisely the convolution integral 

M io-Ps - n7) = - / -— ^ [ImM (O-PS (s) - n7)].ertical cut only (14) 

J (27r)^ 2£;k V ^'4 "^'i n,£.'^J 

provided F (s), the vertex form factor in Im {o-Ps (s) — > 717), is related to the Schrodinger 
momentum wavefunction as (C = \fM jm) 



F (. = 4 (k^ + m^)) = (k) (k^ + m^) = C7<Ao^^ (15) 



Positronium, 




Figure 2 : The vertical cut of the convolution-type factorization. 

Since only the vertical cut is taken into account to get © , this last model obviously cannot 
be analytical in general: Cutkosky rules to get imaginary parts ask for all the possible cuts 
to be included. As we will see, it is precisely a cancellation between vertical and oblique 
cuts that enforces Low's theorem. 

A third way to look at the problem is from the perspective of NRQED scaling rules. 
Typically, the bound state dynamics is soft scale (typical energy < m^- ) while the annihi- 
lation process is hard scale (typical energy > m^^). However, in the corner of phase-space 
where the energy of the emitted photon is very soft, < 7^^ = m? — this separa- 

tion breaks down. Therefore, to get physical energy spectra, it is necessary to keep the 
BE as an arbitrary parameter, and not as a small expansion parameter, during the whole 
computation. 
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2 A new basis for perturbation theory 




Figure 3a : The lowest order loop amplitude. 

From the previous section, it appears that to get analytical amplitudes it is necessary to 
take a four-dimensional loop model (for oblique cuts) and to keep the BE as arbitrary. 
Our proposal for the lowest order amplitude on which perturbation theory is to be built 
is to revert to the Bethe-Salpeter loop amplitude depicted on Fig.Sa, with F (s) given 
in (|15j) . This is a viable lowest order basis because it is easily computed from standard 
point-like QED loop amplitude through differentiation: 



d 




dF 




Figure 3b : Coulomb bound state decay amplitudes from point-like QED ones. 

To see this, it suffices to write the (unsubtracted) dispersion integral for the point-like and 
bound state loop decay amplitudes: 



Point-like case: 



1, (M2) 



1 



ds 



M2 



ImXp (s) 



I W'' ds 

Bound state case: T (M^) = — / ---77 F is) Im Jj, is) 



Thus, with (CSl) 



d 



(16) 



(17) 



Since QED amplitudes are analytical and gauge invariant, so are our lowest order bound 
state decay amplitudes. 

3 Application to QED bound states 
3.1 p-Ps-^ 77 

The point-like amplitude shown on Fig. 4 is written 



Mr 



8m^ 
M2 



e'"''''h,ph,ae*^ ih) 4 ih) Ip with Ip 



-2i 



167r2 



arctan 



'W 



(18) 
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+ Crossed Process 



Figure 4 : The point-like amplitude for p-Ps — > 77 

The corresponding parapositronium decay amplitude is obtained by replacing Zp by I 
obtained from p7j) . The decay rate is then 



r {p-Ps - 77) = arctan - J 



(19) 



Note that this result is non-perturbative with respect to 7. Interestingly, the convergence 
of the perturbative series is greatly accelerated when BE effects are factored out. Writing 

rp-Ps = ^3^(^-arctan-j (l + 5r;.p,) (20) 

we get for ST'^ p^ the same expression © but with much reduced coefficients A'p 0.5326, 
B'p w 0.607 and ^ -3.919. 



3.2 p-Dm— > 7e+e 





Figure 5 : (a) The point-like amplitude for p-Dm — > je'^e , (b) the vertical cut 
and (c) the oblique cut (crossed process understood in each case). 

The simplest process in which the implementation of Low's theorem can be analyzed is the 
annihilation of a paradimuonium (^Sq {fi^fi~) bound state) into a photon and a Dalitz pair 
(see Fig.5a). Following the same steps as above, we get the total rate as the phase-space 
integral over the photon energy spectrum: 



r (p-Dm 



e e 



-7) 



1 — Qe 



1^ dx 


J{x) 




X 



P [X, ae 



(21) 



with 



p (X, fle) W 1 - -7-^ [2 {1~ x)+ fle] -3 

V i — X (1^2^) 

1 /2 



Jix) 1/2 ^ M 47y 

= — — arctan 

X X \TT 27 ttM 

2 \ -1/2 

- 1 



arctan y 



4to^ 



AP (1 



X 



(22) 
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The first (second) term of J {x) jx can be traced to the vertical (oblique) cut, see Fig.5b(5c), 
respectively. 

Low's theorem requirement is met since 



M (p-Dm ^ e+e-7) ° O (u;^) ^ ^ O (1) 

X 



(23) 



To get this behavior, the loop is essential, since taken alone the vertical cut leads to 
J {x) jx ~ The non-perturbative treatment of 7 is equally essential: the limits 

7^0 and x — > are mathematically incompatible 



^^=°^ + i + 0(7) 
X 872 ^ 2 ^ 

This incompatibility is also apparent on the plot of J (x) in Fig. 6. 



(24) 




0.2 0.4 0.6 O.f 



Figure 6 : The behavior of the Coulomb form factor J' (x) , normalized to 
[i7 (3;)]vert cut ' ^ function of X = 2u>^/M for various binding energies 7. 

Finally, note that the total rate is well behaved as 7 — > (as is the surface in Fig. 6). For 
total rate, oblique cuts are subleading. 

3.3 o-Dm^ e+e" 




o-Dm 



Figure 7 : The point-like amplitude for o-Dm e e . 

Another application of interest as a laboratory for leptonic decays of quarkonium is o- 
Dm — >■ e~^e~ as shown on Fig. 7. Proceeding as before, we find 

aM 



r (o-Dm — > e e 



1 + ^ VT^|n(A/2 



3 V- 2,-- 

with n (Af2) obtained from the photon vacuum polarization function Hp (M^) through 
H17|) . Again, we find that the bulk of radiative correction is accounted for at our lowest 
order 



r [o-Dm 



6 



ma 



4- + ... 

vr / NRQED 



(26) 



1 + — VI - fle 1 - 5.33- + .. 
2 / V vr 



Our lowest order 
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This process is especially interesting because the derivative of the photon vacuum 
polarization function obey an anomalous Ward identity [Tj 

A (p') {p'Y - g"'p') = I dxdye'" {a\T {6% (x) J" (y) J- (0)}| 0) 

where 0^ is the trace of the improved energy-momentum tensor. The anomalous term 
— e^/127r^ account for 25% of the BE correction — 5.33a/7r. Obvious open questions are 
whether this is a genuine contribution, and whether perturbative approaches like NRQED 
can catch it. 

3.4 o-Ps— > 777 




The next application is orthopositronium to three photons (Fig. 8a). The computation is 
more involved, but proceeds as previously by taking the derivative of the light-by-light 
scattering amplitude. Details can be found in [21 . The behavior of the photon energy 
spectrum can now be studied as 7 varies (see Fig.Sb): 

* 1 ¥^ ^ '■ ^1 {^^^ + near xi = (correct analytical behavior) 

• 7 = 0: If 7^0 is taken before xi ^ : oc xi (Ore-Powell) 

. 7 ^ 00 : Euler-Heisenberg limit: Ce-h = ^{{F^uFn^ + liF^uF^^n^) 

We observe the same mathematical incompatibility between the two limits uj^ ^ and 
7 ^ as in p-Dm "fe^e~ . However, note well that because of the high correlation in 
the three-photon phase-space, the Low's theorem suppression at low energy modifies the 
high end of the spectrum also. 

For total rate, we write the BE corrections as 

To-Ps = ma'^-^-^B h/M)^^^ (28) 
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with, for small 7, the expansion (note well the slow convergence) 

2 2 s 

^ T T Q. Oi a 

B (7/M) = 1 - 15.4^ + 122^ - 889^ + ... = 1 - 12.1- + 80.2-^ - 502-3 + ... (29) 

and we find again that the bulk of radiative corrections are accounted for as BE corrections: 

2 (tt^ - 9) 

To-Ps = ma^^^—^B (7/M) (l + 6Tlp,,) (30) 

with the same series but with much reduced coefficients A'^ ~ —1.81, B'^ ~ —13.7 and 
C'g ~ 1.48. By the way, note that our method permits an estimation of the yet unknown 
non-logarithmic O (a^) term. 

3.5 Other applications 

Application of the method to spherically symmetric wavefunction is immediate. For exam- 
ple, radial excitation decay amplitudes can be obtained from the punctual QED amplitude 
as 



n 



iFi ( 1-n, 2,167^ ^ 



9M2 



7n (M^) (M2) 



(31) 



where the hypergeometric function is essentially the Laguerre polynomial. 




Figure 9 : The point-like mass renormalization graphs. 

Hyperfine splitting can also be dealt with using the present method. From double deriva- 
tives of the punctual mass renormalization amplitudes depicted on Fig. 9 

1 \ / d 

^para{ortho) {M"^) = {32ttC (l)o7) 2 \'dM^ ) ^P'P"-Mortho) {M'^) (32) 

and taking into account the one-photon annihilation graph, we find again that our low- 
est order AEhf = {0.53ma'^) lowest order reproduces much of the NRQED correction^ 
AEhf = (0.5833ma^)^j^Q^^. For further details, see 

4 Application to QCD bound states 

The charmonium and bottomonium are the strong analogues of the positronium. Their 
binding energies are relatively small 

Positronium (e"e+) Aml/M^ wl + 10"^ (ss 1 + 0^/4) 
Charmonium (cc) 4m|,/Mj^^ 1.45 (33) 
Bottomonium (66) Am\/My 1.25 
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Their decay rates are computed using adapted to the present case 



2S+lg 

and we have the analogies (JZ 



^) « \4>o\'^ {vrelCfscatt. {QQ ^ X)) 



/M^ and Cc the electric charge of the c) : 



Table I : Correspondence between QED and QCD applications. 



QED 



QCD 



p-Ps 77 



o-Dm 
o-Ps - 



-> e^e 
777 



Vc 77 

Vc gg 

J/^ ^ 7* - 
J ^ 777 
J/V' ^ 7ra 

Jl"^ ^ ggg 



e ' e 



Decay ratei^S; 



f (7r2 - 9) e6a37^ 
^ (tt^ - 9) e2aa|7^ 

# (vr^ - 9) 4^ 



(34) 



We now review a number of observables in quarkonium physics that may be affected by 
BE effects. Since we do not know the precise form of the quarkonium wavef unctions, the 
following discussions are rather qualitative. 



4.1 Prompt photon spectra 

The vector quarkonium differential decay rates into three gauge bosons are: 



64 K o 
777) = y ^ (^) 



dx 

dr 

dx 



igg) 
ggg) 



128 



160 



enaogTZ (x) 



ai^n n (x) 



(35) 



with the Ore-Powell spectrum function (x) given in (|11|) . which contradicts Low's theo- 
rem requirements. In Sec. 3. 4 it was shown that BE effects introduce a softening of the spec- 
trum at high energies, and this is precisely what is observed for both J/V' ^ 7 + hadrons 
and T ^ 7 -|- hadrons, compare Fig. 8b and Fig. 10 (where the low energy increase is due 
to bremsstrahlung processes ^J). 



1500 



N.„.dx 



■ Background Subtracted Data 
I — I Lowest order QCD (Ore-Powell) 



■ Background Subtracted Data 
Lowest order QCD (Ore-Powell) 



1000 




dN 
dx 



500 




0.25 



0.50 



0.75 



1.00 



Figure 10: Photon spectrum in J/-^ ^ 7 -|- had. 9^ and in T ^ 7 -|- had.\\()\. 
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4.2 pyr puzzle 

The decay rates in table I are valid for both J/ifj and i/^(25'), the only change being in 
the wavefunction (pQ and mass M. From them, one could expect that the ratio of the 
three-gluon decay rates of J/V' and -0 (25) is the same as the one into e+e". Going one 
step further, one could argue that the ratio into exclusive hadronic modes would also be 
constant ^ : 

B{^{2S)^X) ^ B{i^{2S)^e+e-) 
B{J/iP^X) ^ B{J/ij^e+e-) ^ ' 

Experimentally, the ratio of leptonic modes is roughly 14%, while for hadronic modes, 
some are close to 14% and others much suppressed 13 : 

The BE corrections to the one-virtual photon and three-gauge boson rates originate 
in the dynamics of the photon 2-point and 4-point functions, respectively. Those two 
are completely different: the 4-point function is much more sensitive to the ratio of the 
masses of the loop particle and initial particle. As a result, no matter the form of the 
wavefunction, three-gauge boson modes will be much more affected for non-negligible 7. 
Now, the isospin violating mode cjtt proceeds through a virtual photon, as is the leptonic 
mode, and thus both receive the same correction. On the other hand, the 14% rule will 
be invalidated by large BE correction for the pir mode, which proceeds from ggg. 

4.3 Extraction of as and the perturbation series 

The strong coupling constant can be extracted from ratios of inclusive hadronic modes 
and electromagnetic ones|13|. For instance, from r/^ and J/ip modes, one gets 

r]c 77 

Vc gg 




as (me) « 0.31 Jli)^^gg \ ^Qs(me)« 0.19 (38) 

as(Mz)«0.12 J/^^ggg ] ^as(Mz)«0.10 



For r/c, both decay modes involve the same dynamics and their ratio is unaffected by BE 
corrections; as is in agreement with the world average |T3j as [Mz) ~ 0.117 it 0.002. On 
the other hand, for J /ip decays, the dynamics are very different and large BE corrections 
may arise. 

As we have repeatedly seen in the case of positronium, the possibility of factoring out 
BE corrections improves the behavior of perturbative series. In the context of quarkonium, 
where first order corrections tend to be of the same order as the leading term, this may 
be crucial to get meaningful results. To take the worst casejH] 

2_.^ ,6 



r ( J/V' ^ 777) = Y K - 9) (1 - 12-6^ j (39) 

For as ~ 0.25, the width vanishes. Now, the BE correction can be sizeable 

r ( j/^ ^ 777) = y - 9) B [ij/^/Mj,^) (40) 

If the charmonium wavefunctions were Coulombic, B {jji^jMjj^ « 0.01, illustrating 
that really large suppression may indeed occur. 



10 



5 Conclusions 



For QED bound states, a new basis for perturbation theory is proposed, in agreement with 
analyticity. Differential rates then behave as predicted by Low's theorem. Also, BE effects 
can be factored out of the perturbation series, thereby accelerating their convergence. 

For QCD bound states, a perturbative expansion in 7 like in NRQCD appears as 
very questionable, see the slow convergence of H19|) or ()29() . BE effects have to be dealt 
with non-perturbatively, and can change the whole picture. In particular, prompt photon 
spectra in vector quarkonium decays are softened. Also, the difference in the dynamics of 
the photon 2-point and 4-point function can explain, at least in part, both the /ovr puzzle 
and the smallness oi as as extracted from vector quarkonium. 

Work is in progress to extend the method to higher order computations, P-wave 
positronium and quarkonium decays, transitions among states, production rates,... 
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